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Abstract 

Cool objects glow in the infrared. The gas and solid-state species that escape the 
stellar gravitational attraction of evolved late-type stars in the form of a stellar 
wind are cool, with temperatures typically ^1500 K, and can be ideally studied 
in the infrared. These stellar winds create huge extended circumstellar envelopes 
with extents approaching 10^^ cm. In these envelopes, a complex kinematical, 
thermodynamical and chemical interplay determines the global and local struc- 
tural parameters. Unraveling the wind acceleration mechanisms and deriving the 
complicated structure of the envelopes is important to understand the late stages 
of evolution of ~97% of stars in galaxies as our own MiUcy Way. That way, we 
can also assess the significant chemical enrichment of the interstellar medium by 
the mass loss of these evolved stars. The Herschel Space Observatory is uniquely 
placed to study evolved stars thanks to the excellent capabilities of the three in- 
frared and sub-millimeter instruments on board: PACS, SPIRE and HIFI. In this 
review, I give an overview of a few important results obtained during the first two 
years of Herschel observations in the field of evolved low and intermediate mass 
stars, and I will show how the Herschel observations can solve some historical 
questions on these late stages of stellar evolution, but also add some new ones. 
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1. Introduction 

The launch of the Herschel Space Observatory (Pilbratt et al. 2010) in May 
2009 can be marked as a milestone in the area of infrared astronomy. With a 
mirror of 3.5 m in diameter, Herschel studies the Universe in the far-infrared and 
sub-millimeter wavelength ranges (60-670 yum) thanks to the capabilities of the 
three instruments on board: HIFI (de Graauw et al. 2010), PACS (Poglitsch et al. 
2010), and SPIRE (Griffin et al. 2010). The infrared (IR) and sub-millimeter 
regions of the spectrum are of great scientific interest, not only because it is here 
that cool objects (10-lOOOK) radiate the bulk of their energy, but also because 
of its rich variety of diagnostic atomic, ionic, molecular and solid-state spectral 
features. Measurements at these wavelengths permit determination and evaluation 
of the physical processes taking place in astronomical sources, estabUshing the 
energy balance, temperature, abundances, density and velocity. 

The Herschel instruments enable the astronomers to study the coldest and 
dustiest objects in space, among which star forming regions, young stars, evolved 
stars, dusty galaxies, ... In this review, we will focus on the late stages of evolu- 
tion of low and intermediate-mass stars (LIMS, M~0.6-8Mo), the mass interval 
to which most of the stars in the Universe belong. Using a standard Salpeter initial 
mass function (IMF), 97% of all stars bom will go through the so-called Asymp- 
totic Giant Branch (AGB) phase, where mass loss dominates the stellar evolution. 
More massive stars will lose a significant fraction of their mass in the so-called 
supergiant phase. The mass loss of AGB stars dominates the total mass of the 
interstellar medium (ISM) and almost half of the heavy elements returned to the 
ISM originates in the stellar interiors of these stars. 

The detailed physical description of the mass-loss process for late-type stars is, 
however, still unknown after several decades of studies. The generally accepted 
idea on the mass-loss mechanism for asymptotic giant branch (AGB) stars is based 
on pulsations and radiation pressure on newly formed dust grains. However, 
oxygen-rich AGB stars suffer from the so-called 'acceleration deficit' dilemma, 
which states that mass-loss rates due to the formation of siHcate dust alone are 
orders of magnitude smaller than observed ones (Woitke 2006). The formation of 
large grains in the inner wind region might solve this riddle (Hofner 2008; Norris 
et al. 2012). The proposed models for the mass loss of AGB stars are very unlikely 
to be applicable to red supergiants (RSGs) since they are irregular variables with 
small ampHtudes. Processes Unked to convection, chromospheric activity or rota- 
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tion might play an important role as trigger for the mass loss in RSGs. Josselin 

and Plez (2007) proposed that turbulent pressure generated by convective motions, 
combined with radiative pressure on molecular lines, might initiate mass loss, but 
also Alfven waves generated by a magnetic field might contribute (Hartmann and 
Avrett 1984). 

Several guaranteed time and open time Herschel proposals (will) focus on dif- 
ferent aspects of these late stages of stellar evolution dominated by mass loss. 
E.g. large maps made with the PACS and SPIRE photometers will probe the in- 
frared signature of dust emitted by evolved stars and supernovae in the Magel- 
lanic Clouds (Meixner et al. 2010). The goal is to assess the life cycle of matter 
in (external) galaxies by deriving general properties of the circumstellar material 
(CSM) surrounding these evolved stars. However, the details on the onset of the 
mass loss, on the mass-loss history, and on the chemistry, kinematics and thermo- 
dynamics in the complex outflows (see Fig. 1) can only be assessed for the most 
nearby evolved stars (closer than few kpc). Unravelling these physical and chem- 
ical details will enhance our knowledge on the mass loss process in evolved stars, 
and will hopefully enable us to derive some general quantities/laws applicable to 
more distant stars. 

In this review, I will give an overview of some of the new scientific insight 
gained thanks to Herschel after being 2.5 yr into mission. The focus of Sect. 2 is 
on new geometrical and dynamical insights in the evolved stars' circumstellar en- 
velopes (CSEs) as created by their stellar winds. Most of the constraints will come 
from PACS and SPIRE images, with photometer bands at 70, 100, and 160//m, 
and 250, 350, and 500 ^um respectively. Sect. 3 illustrates some spectroscopic re- 
sults obtained with the high-spectral resolution capabilities of HIFI (157-212//m 
and 240-625 fxra, R as high as 10^) and the medium-spectral resolution PACS 
(60-210yum, R = 1000 - 5000) and SPIRE (194-672 //m, R = 40 - 1000) spec- 
trometers. In each of these two sections, the examples will cover subsequently 
the different geometrical regions displayed in Fig. 1. In Sect. 4, I will discuss 
the answers on some 'old' astronomical questions concerning these late stages of 
evolution of LIMS, but will also add some new questions to which the answer is 
still unclear. Some future perspectives are given in Sect. 5, and some conclusions 
are dravm in Sect. 6. 

2. New geometrical and dynamical insights from Herschel 

Despite the importance of the mass-loss process in terminating the AGB phase 
and in replenishing the ISM with newly-produced elements, the nature of the mass 
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Figure 1: Schematic drawing (not to scale) of an AGB star. Pulsations and dust formation are the 
two key ingredients for the generation of a stellar wind, with typical terminal velocities between 
10-30 km/s and mass-loss rates between 10"*^ and 10"'*Mo/yr Several chemical processes are 
indicated (in red) at the typical temperature and radial distance from the star where they occur 
['TE' stands for thermal equilibrium]. The penetration of interstellar ultraviolet photons in the 
outer envelope is shown in blue. Understanding the chemical processes in AGBs is crucial to 
know the chemical content of the (unprocessed) ISM. The results described in Sect. 2 and Sect. 3 
will cover subsequently the different geometrical regions displayed in this figure. 
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loss is still not well understood; in particular, one cannot predict the mass-loss rate 
for a star of given properties. 

While for modeling purposes, the envelope is often assumed to be spherically 
symmetric and formed by a constant mass-loss rate, observations show more and 
more evidence contradicting these two basic assumptions. CO and scattered light 
observations of (post) AGB objects show that winds from late-type stars are far 
from being smooth. The shell structures found around, e.g., AFGL 2688 (Sahai 
et al. 1998) and IRC +10216 (Mauron and Huggins 1999) indicate that the outflow 
has quasi-periodic oscillations, with density contrast corresponding to a factor 
~100-1000 (Schoier et al. 2005; Decin et al. 2006). These mass-loss modulations 
on a time scale of a few hundred years may be due to a complex feedback between 
hydrodynamics, gas/dust drift, and gas-solid chemistry (Simis et al. 2001). At the 
tip of the AGB, the mass ejection process becomes very strong, with mass-loss 
rates as high as 10^"^ Mg/yr. 

Herschel data are now adding more and more complexity to this geometrical 
picture. In this next subsections, recent Herschel results will be described cover- 
ing gradually the envelope from the inner to the outermost structures. 

2.1. Inner envelope structure 

The PACS and SPIRE image capabilities for nearby evolved stars are tremen- 
dous. Thanks to their high spatial resolution (of ~V'), one is able to map the 
structure of the inner envelope surrounding these evolved stars in quite some de- 
tail. Deducing if the envelope shows deviations from a smooth ID envelope struc- 
ture, created by a constant mass-loss rate, is however not so straightforward due 
to the complex point-spread function (PSF) of the PACS photometer and the fact 
that the cool central target (with typical temperatures between 2000-3500 K) can 
be quite bright at these infrared wavelengths. 

To emphasize the shell morphology in the inner envelope, one can remove 
the extended envelope halo by subtracting a smooth, azimuthally averaged profile 
represented by a power law r~" (see Fig. 2, Decin et al. 2011). The method has 
been tested in detail (Royer et al., in prep.). An illustration for a point-like source 
is shown in the upper panels of Fig. 3: for the oxygen-rich semi-regular pulsating 
AGB star R Dor, experiencing a low mass-loss rate of ~ 1 .2 x 10"^ Mo/yr (Schoier 
et al. 2004) no indications can be found for deviations from a smooth spherically 
symmetric envelope at scales >15". 

The images of the three other targets shown in Fig. 3 and the image of IRC -i- 1 02 1 6 
(see Fig. 2) testify to a much more complex envelope structure. 
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Figure 2: PACS images of IRC +10216 at 70//m (left), 100 //m (middle) and 160 yum (right) after 
subtraction of the smooth halo of the circumstellar envelope (Decin et al. 2011). The six radial 
spikes are due to the support spider structure of the secondary mirror The field-of-view is 16'xir. 

IRC +10216. is the nearest carbon-rich AGB star at a distance of ~ 150 pc (Crosas 
and Menten 1997). Mauron and Huggins (1999) have detected multiple, almost 
concentric, shells (or arcs) out to a distance of 55" in radius on top of the smooth 
extended envelope. The Herschel images show the discovery of multiple dust 
shells until 320". Decin et al. (201 1) argue that the origin of the shells is related 
to non-isotropic mass-loss events and clumpy dust formation. 

X Cyg. is the prototype of a S-type (C/0 ratio slightly lower than 1) Mira vari- 
able with a period of 408 days at a distance of ~150 pc and a mass-loss rate of 
~ 5 X 10~^Mo/yr (Schoier et al. 2011). The PACS image shows some clear in- 
tensity enhancements in the north-western and south-eastern direction. Ragland 
et al. (2006) reported the detection of asymmetric brightness distributions using 
the IOTA interferometer in the //-band, with a angular resolution of 5-10 mas, 
probably arising from just above the stellar photosphere. This is confirmed by 
Tatebe et al. (2006), who interpreted their interferometric data in terms of a steady 
excess of dust emitted to the east, or alternatively, as being due to a hot spot on 
the east side of the star that illuminates a symmetric dust shell in an asymmetric 
way. According to Tatebe et al. (2006) a companion is unlikely to be the cause 
of the reported asymmetry. The Herschel/HIFI line profiles of the H2O, SiO and 
HCN show a consistent weak asymmetry, which are attributed to small-scale devi- 
ations from spherical symmetry and the presence of a (weakly) clumped medium 
(Schoier et al. 2011). 

fi Cep. is an oxygen-rich supergiant at a distance of 390 pc (Cox et al. 2012). The 
PACS 70 jum image is dominated by the bow shock structure situated at ~50" (see 
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Figure 3: PACS images of four evolved stars before (left) and after (right) subtraction of the 
smooth halo of the circumstellar envelope. North is to the top, east is to the left. The six radial 
spikes (visible in the image of R Dor) are due to the support spider structure of the secondary 
mirror The red circle with a radius of 15" marks the region where some PSF artifacts are still 
visible. These data were obtained in the framework of the Herschel MESS Guaranteed Time Key 
Programme (Groenewegen et al. 2011). A full analysis of the inner envelope structure of some 
tens of stars will be presented by Royer et al. (inJprep.) 



Sect. 2.3). Some asymmetries at ~20" away in the envelope are seen in the PACS 
image, although it is not clear if this is related to a turbulent mass-loss history or 
due to the large instabilities in the bow shock. 

a Ori. is also an oxygen-rich supergiant at a distance of 197pc (Harper et al. 
2008). Many detailed studies have already described some characteristics of the 
complex atmosphere, chromosphere and dusty envelope of Betelgeuse. Irregular 
structures have been reported within 2.5" from the central target (e.g. Kervella 
et al. 2011). The Herschel images show the first evidence for a high degree of 
dumpiness of the material lost by the star beyond 15", which even persists until 
the material collides with the ISM (Decin et al. 2012a). Very pronounced asym- 
metries are visible within 110" from the central star, although some weaker flux 
enhancements are visible until ~300" away. The typical angular extent ranges 
from -10-90°. 

Indeed, five examples of evolved AGB and supergiant stars, and five differ- 
ent mass-loss histories. There is one common denominator: the emission we see 
is (almost completely) caused by dust - there might be some minor contribution 
from some atomic or molecular lines in the photometer bands as [O I] at 63 //m. 
This means that we are the direct witnesses of the role of dust in creating the enve- 
lope structure around these evolved stars. In the Herschel/MESS GTKP some 78 
AGB and supergiant stars are observed; other open and guaranteed time proposals 
will image some dozens of other nearby evolved stars. It is clear that these results 
will revolutionize our ideas on the wind creation in these targets, and that these 
images form the perfect ground for follow-up observations with ALMA. 

2.2. Wind acceleration 

The formation of dust species is thought to play a crucial role for the onset 
of the wind acceleration in AGB stars. In case of the supergiants, the role of the 
dust species for the wind acceleration is absolutely not clear. But even in case 
of AGB stars, it is not yet understood how the wind can be driven in an oxygen- 
rich environment, where the glassy character of the oxides and pure (magnesium- 
rich) silicates — the species which might form close to the star — prevents them 
from gaining enough momentum (Woitke 2006) ^ The formation of both carbon 
and silicate grains (Hofner and Andersen 2007) or micron-sized Fe-free silicates 



When solving the momentum equation shown in Fig. 4, the silicate dust composition as de- 
rived by Justtanont and Tielens (1992) was assumed. The main component is amorphous oUvine 
(Mg;,Fei_;,Si04). 



8 



(Hofner 2008; Norris et al. 2012) have been proposed as possible solutions to this 

dilemma. 

The high-spectral resolution of HlFl enables us to determine the expansion ve- 
locity of the spectral lines which sample different excitation regimes and hence are 
formed in different regions in the envelope. In addition, spectral lines (partially) 
formed in the wind acceleration region have a typical Gaussian profile (instead of 
a flat-topped, two-horn or parabolic profile, Lamers and Cassinelli 1999). Dedi- 
cated modelling can then yield some observational evidence on the acceleration. 

In the classical jS-parametrization (e.g., Lamers and Cassinelli 1999) used to 
approximate the velocity structure. 



with vo the velocity at the dust condensation radius and Uoo the terminal velocity, /3 
quantifies the wind acceleration. Schoier et al. (201 1) derived a value ofp = 0.95 
for the S-type AGB star;^f Cyg; for the oxygen-rich AGB star IK Tau a value of 1 is 
derived (see Fig. 4 from Decin et al. 2010d). These are significantly more shallow 
velocity laws compared to that derived by Schoier et al. (2006) for the high-mass- 
loss rate carbon-rich AGB star IRC -1-10216, where yS = 0.2. This illustrates that 
the wind driving and momentum coupling is indeed more efficient in carbon-rich 
AGB stars, but also that oxygen-rich and S-type AGB stars do reach their terminal 
velocity within ~100Rik.. 

An intriguing figure on observed expansion velocities in oxygen-rich AGB 
stars is presented by Justtanont et al. (2012) (see Fig. 5). In that figure, a function 
similar to Eq. 1 has been overplotted, but with a dependence on the energy rather 
than the radius: 



where Eo is an arbitrary initial energy. Five targets show a clear wind acceleration: 
TX Cam (J3 = 0.2), R Cas (J3 = 0.5), IK Tau (fi = 0.9), o Cet (J3 = 0.9), and 
IRC +10011 (jS = 0.9); in W Hya and AFGL 5379 an increasing velocity w.r.t. 
the upper energy level is visible, although a strong onset of the acceleration is 
absent. In the case of R Dor, however, the highly excited lines of both H2O and 
SiO indicate a larger expansion velocity in the inner part of the envelope than for 
the lower J-transitions of CO. 

It is important to realize that the interpretation of this kind of figure is not 
unequivocal since different molecules are involved, each having different optical 
depths and different Une formation regions, and some lines are even masing. A 
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Figure 4: Velocity profile of IK Tau as published in Decin et al. (2010d). Velocity data are obtained 
from mapping of maser emission: SiO (Boboltz and Diamond 2005), H2O (Bains et al. 2003), and 
OH (Bowers et al. 1989). The CO expansion velocity derived from ground-based CO J=l-0 data 
is also indicated (Decin et al. 2010c). The triangles show the place in the envelope where the 
line formation is highest for the HIFI data presented in Decin et al. (2010d), including '^CO, 
i^CO, H'^-O, H'^O, Hi'^O, 28SiO, ^^Si), ^°SiO, HCN, and SO lines. The horizontal bai's show the 
minimum and maximum radial distance for the line formation of each individual transition. The 
vertical bars show the uncertainty on the observed line widths. The expansion velocity deduced 
from solving the momentum equation (assuming extinction efficiencies for spherical dust particles 
with a main component being amorphous olivine) is shown by the full green line. The dashed 
blue line represents a power law (Eq. 1) with /? = 1. For comparison, an even smoother expansion 
velocity structure with y6 = 1.8 is shown with the red dotted line. The vertical dashed black line 
indicates the dust condensation radius Rinner- The velocity at Rinner is assumed to be equal to the 
local sound velocity. 
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Figure 5: Observed expansion velocity for different molecular lines observed with HerschelfHWl 
in a sample of oxygen-rich AGB stars (Justtanont et al. 2012). The abscissa gives the upper energy 
level of the different lines, the ordinate the measured expansion velocity. 
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varying turbulent velocity might play a role to explain the velocity profiles show- 
ing an unexpected deceleration w.r.t. the upper energy level, but cannot explain 
the total spread. In general, the turbulent velocity in the outer atmosphere and 
inner envelopes of AGB stars is not well known. For non-Mira M-giants, the 
analysis of CO first-overtone bands indicates a turbulent velocity ranging from ~1 
to 4 km/s (Tsuji 1986). For the well-known carbon-rich AGB star IRC +10216, 
it was speculated that the turbulent velocity decreases somewhat with increasing 
distance from the star, from a value of a few km/s close to the stellar surface to 
~0.6 km/s at hundreds of stellar radii (Skinner et al. 1999). Using interferometric 
observations, a turbulent velocity of (only) 1.5 km/s in the inner wind of R Dor 
was derived (Schoier et al. 2004), far less than the spread seen in Fig. 5. But also 
the high-resolution HIFI data offer the possibility to derive the turbulent velocity 
from a detailed analysis of the line wings. A recent example showing the strength 
of the HIFI data to estimate the turbulent velocity is shown in Fig. 8 in De Beck 
et al. (2012). This strength of the HIFI data is still largely unexplored. 

It is clear that from a detailed analysis of the high-spectral resolution HIFI data, 
one can derive the global (ID) characteristics of the velocity field after a dedicated 
analysis of different lines with different excitation levels. However, line formation 
regions are usually quite broad (see Fig. 4), complicating a precise determination 
of the wind acceleration. 

2.3. CSM-ISM bow shock 

The Herschel PACS and SPIRE images provide abundant scope for detailed 
studies on the complex interaction region between the circumstellar and interstel- 
lar material. 

Photometric observations performed with the IRAS telescope (Noriega-Crespo 
et al. 1997), and later confirmed with the Spitzer Space Telescope and AKARI 
(Martin et al. 2007; Ueta et al. 2006, 2008) indicated a new type of process that 
might impact the chemical composition of circumstellar material (CSM) which 
is injected into the ISM: the interaction zone between the CSM and ISM might 
be more energetic than assumed earlier, resulting in a bow shock structure at the 
place where the wind collides with the ISM (in analogy with the solar environ- 
ment). However, these early observations of the far-IR and UV-bright emission 
structures lack any spectral resolution, and as they have only poor spatial resolu- 
tion, allow only basic morphological studies (e.g. Wilkin fitting (Wilkin 1996) to 
derive relative velocities). 

In the MESS GTKP (Groenewegen et al. 201 1) several known and new interac- 
tion objects have been imaged with the PACS instrument, revealing unprecedented 
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Figure 6: PACS 70/im images of the CSM-ISM interaction of a sub-sample of evolved stars as 
observed in the MESS GTKP. More details and the image of other targets can be found in Cox 
et al. (2012). The images are best viewed on screen. 



detail. A nice overview is given by Cox et al. (2012) (see also Fig. 6). The detec- 
tion rate of bow shocks is very high (80%) when the subset is hmited to distances 
less than 300 pc! The detached shells around TT Cyg, U Ant, and AQ And are 
discussed by Kerschbaum et al. (2010) and on some more objects by Kerschbaum 
et al. (201 1). Bow shocks are reported for CW Leo (Ladjal et al. 2010), X Her and 
TX Psc (Jorissen et al. 201 1), and a Ori (Decin et al. 2012a). For o Cet the inner 
part of the stellar wind bubble bounded and formed by the termination shock is 
seen (Mayer et al. 2011). 

The morphological classification introduced by Cox et al. (2012) includes 'fer- 
mata's' (cfr. X Pav), 'rings' (cfr. U Hya), 'eyes' (cfr. VY UMa), and 'irregular'- 
like structures. Depending on the space velocity, mass-loss rate, ISM temperature 
and density, etc. different bow shock structures might appear (see Fig. 8-10 in Cox 
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Figure 7: PACS 70/im image of R Leo (left) compared to a hydrodynamical simulation of the 
wind-ISM interaction computed with the AMRVAC code (van Marie et al. 2011). The parameters 
for this simulation are a stellar wind velocity of 15km/s, a constant gas mass-loss rate of 1 x 
lO^^Mo/yr, a dust-to-gas mass ratio of 0.01, a space velocity of 25km/s, a local ISM density 
of 2 cm"^, and an ISM temperature of 3 K. The upper right figure shows the gas density, which 
range (in log-scale) between 10"^'* and 10 g/cm^; the lower right figure represents the dust grain 
particle density, ranging (in log-scale) between 10"'" and lO^^*^ cm"^. 
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et al. 2012), with clear signatures of Kelvin-Hehnholtz^ and/or Rayleigh-Taylor^ 

instabilities. A textbook example in this respect is R Leo (see Fig. 7), where clear 
'smaller scale' Rayleigh-Taylor instabilities (visible in the form of mushroom- 
like signatures) and 'large scale' Kelvin- Helmholtz (deforming the RT instabil- 
ities along the contact discontinuity) are present. These Herschel images are the 
first images showing both types of instabilities in the shocked wind-region around 
evolved AGB and supergiant stars. 

The new Herschel data trigger the development of hydrodynamical codes to 
simulate these environments (e.g. van Marie et al. 2011; Mohamed et al. 2012). 
Fig. 8 shows the wind-ISM interaction region obtained from a hydrodynamical 
simulation using the AMRVAC code (Keppens et al. 2012). This code has recently 
been improved (van Marie et al. 201 1) to include dust grains in the stellar wind, 
and to take the drag forces between dust and gas into account. The example shows 
that the dust is only heated up by a few K in the interaction zone, while the gas 
temperature can be as high as 10 000 K. Clumps in the inner envelope (as seen, 
e.g., in the Herschel images described in Sect. 2.1) might have an imprint on the 
bow shock morphology. Simulations show that a smoothly varying (sinusoidal) 
wind density, as shown in Fig. 9, has an imprint in the free flowing wind, but 
that the density contrast is completed erased at the turbulent contact discontinuity. 
Only very strong 'picket-fence '-like variations with a factor of a few thousand 
might impact the bow shock morphology (Decin et al. 2012a). 

The 'ring' -like bow shocks (as U Hya) all belong to thermally-pulsating stars, 
being carbon-rich AGB stars or technetium enriched M-type supergiants. The 
origin of the detached shells appears linked to the interaction of a fast wind with 
an 'older' slow wind or a dramatic change in mass loss caused by, e.g., a thermal 
pulse (TP). 

The structures least understood, in my opinion, are the 'eye'-like structures 
(as VY UMa), two elliptical non-concentric arcs observed at opposing sides of 
the central source with a covering angle of <180°. For a few sources, the two 
arcs are connected and there is even tentative evidence for a jet structure in the 
mid plane. Binarity might play a role, but this can not be confirmed or excluded 
at the moment. Magnetic fields might play a role as well, but this field is still 
totally unexplored. As an example, we show in Fig. 10 a simulation using the 



^Kelvin-Helmholtz (KH) instabilities arise when velocity shear is present in a fluid or when 
there is sufiicient velocity difi'erence across the interface between two fluids 

^Rayleigh-Taylor (RT) instabilities might occur at the interface between two fluids of diff'erent 
densities, when the hghter fluid is pushing the heavier fluid 
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Figure 8: Hydrodynamical simulation of the wind-ISM interaction of the envelope around Betel- 
geuse after 50 000 yr of evolution using the AMRVAC code (van Marie et al. 20 11 ). The parameters 
for this simulation are a stellar temperature of 3650 K, a stellar wind velocity of 15km/s, a con- 
stant gas mass-loss rate of 3 x 10"^ Mo/yr, a dust-to-gas mass ratio of 2.5 x lO"-', a velocity w.r.t. 
the local ISM of 28.3 km/s, a local ISM density of 3.2 x lO*"^"*' g/cm^ and an ISM temperature 
of 3 K. First panel: gas density in glcvc?, second panel: dust density in gjcvp? (for a dust parti- 
cle with a radius of 0.05 yum), third panel: gas temperature in Kelvin, fourth panel: gas absolute 
velocity in cm/s. After a simulation time of 50000yr, the place of the bowshock interaction has 
stabilized: the termination shock occurs at ~0.33pc, the bowshock at ~0.46pc, and the turbulent 
astropause/astrosheath in between. The instabilities are Rayleigh-Taylor instabilities (due to the 
density difference between the CSM and ISM material) which are somewhat 'dis-formed' due to 
the velocity shear between CSM and ISM (yielding Kelvin-Helmholtz instabilities). The gas tem- 
perature in the interaction zone can be as high as few thousand K, while the dust temperature only 
increases slightly w.r.t. the region without ISM interaction. Detailed inspection of the gas and dust 
density shows that a small dust particle (of 0.05 //m in this simulation) is able to travel beyond 
the CSM-ISM contact discontinuity, but does not enter the unshocked ISM. A dust particle with a 
grain size of 1 jum is able to cross the border between shocked and unshocked ISM material. 



16 



Figure 9: Same hydrodynamical simulation as in Fig. 8, but for a varying mass-loss rate with a 
period of 1000 yr, implemented as M(f)= 1.8 x 10"^ + (1.5 x 10"^ x sin(2;rf/1000yr))MQ/yr. 

MPI-AMRVAC code (Keppens et al. 2012) for a star, on which a toroidal magnetic 
field 



(with Vrot the surface rotation of the star, i?* the stellar radius, r the radial grid 
coordinate, and 6 the colatitudinal angle) has been imposed. The stellar magnetic 
field 5* is defined as a fraction of the kinetic energy of the wind (Begelman and 
Li 1992): 



MVoo \Vool 

This setup for the magnetic field is similar to the one used by Garcia-Segura et al. 
(1999), except that Eq. 3 reverses the polarity of the magnetic field at the equator. 
The stellar parameters are 7?* = ISORq and Vrot = lOkm/s. The influence of rota- 
tion on the stellar wind was not included. The simulation was run for cr = (no 
magnetic field) and cr = 0.1 (yielding a stellar magnetic field B^, = 0.931 Gauss)"*. 



The detected longitudinal field in the red supergiant Betelgeuse is ~1 Gauss (Auriere et al. 
2010) 
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This simulation shows indeed an 'eye'-like shape, but we should take into account 
that no magnetic fields are yet detected in AGB stars. 
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Figure 10: Logarithm of the density in g cm"^ for the expansion of an AGB wind into a constant 
density ISM; one without magnetic field (cr = 0, /e/f) and one with B - 0.93 1 G (for cr = 0. 1, right). 
The non-magnetic simulation shows a very thin shocked wind layer and a thicker shocked ISM 
shell. For the simulation with magnetic field, both layers are approximately equally thick. The 
magnetic field also creates collimated streams of matter over the poles and at the equator, which 
distort the spherical symmetry of the shell. 

I am convinced that these Herschel images are opening a new era of detailed 
studies of this intriguing interaction zone between ISM and CSM. Highly sophis- 
ticated model simulations and dedicated observations (with, e.g., ALMA) will 
elucidate the key astrophysical quantities triggering the appearance of this whole 
zoo of interaction phases. 



3. New chemical insights from Herschel 

To date, more than 70 molecules have been detected in the circumstellar en- 
velopes (CSEs) of evolved stars (e.g. Cemicharo et al. 2000; He et al. 2008). 
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Figure 1 1 : Comparison between the rich HerschelfPACS spectra of the O-rich AGB star R Dor 
(black), the O-rich supergiant VY CMa (gi-ey) and the C-rich AGB stai- CIT 6 (red). Spectra of 
VY CMa and CIT 6 are scaled to facilitate the comparison. Most lines are blends, with he most 
prominent features being CO and H2O transitions (O-rich) and HCN (C-rich). Each spectrum 
contains some few hundred lines. Line identifications can be found in Decin et al. (2010a,b); 
Royer et al. (2010). 



Admittedly, most of them are only detected in the nearest carbon-rich AGB star 
IRC +10216. The effect of nucleosynthesis in the stellar core and subsequent con- 
vective envelope mixing or dredge- up s determine the abundance structure of the 
stellar outer atmospheric layers. Thermodynamic equilibrium and non-equilibrium 
reactions, photochemical reactions, ion-molecule reactions, and the condensation 
of dust grains establish the abundance stratifications throughout the circumstel- 
lar envelope (see Fig. 1). The IRAS, Spitzer and AKARI data, and recently the 
Herschel PACS and SPIRE images (see Sect. 2.3) indicate that the CSM-ISM in- 
teraction region might be a chemically complex region. 

The three spectrometers on board Herschel have proven to be three powerful 
instruments to unravel these complex chemical characteristics. Some first results 
are given in the subsequent sub-sections covering gradually the full envelope start- 
ing at the innermost regions. It will become clear that most of these first results 
are based on thorough non-local thermodynamic equilibrium (non-LTE) calcula- 
tions, sometimes invoking the modeling of different molecules at once. Indeed, in 
the era before Herschel, often only a handful of lines (at maximum) of one spe- 
cific molecule were modeled by solving the statistical equilibrium rate equations. 
Today, often hundred molecular line transitions of several molecules are simulta- 
neously modeled (see Fig. 11), if needed by using supercomputers. In addition, it 
becomes more and more clear that a consistent dust-gas modeling effort is of ut- 
most importance to correctly derive the physical and chemical envelope structure 
(see, e. g., Lombaert et al. 2012). 

3. 1. Nucleosynthesis and dredge-ups 

The material injected in the ISM reflects the effects of nucleosynthesis in the 
stellar core and subsequent convective envelope mixing or dredge-ups (DUs). The 
effect of nucleosynthesis and DUs on the stellar abundances can be constrained by 
different observations: one can study the abundance pattern (1) in the dense outer 
layers of the stellar AGB atmosphere, (2) in the tenuous extended CSE, and (3) 
from presolar grains which are thought to be condensed from AGB stars. How- 
ever, stellar atmosphere analysis excludes high mass-loss AGB stars and presolar 
grains have unpredictable selection effects. Due to these observational selection 
effects, the injected abundance pattern is best studied in the stellar wind of the 
AGB or supergiant star. 

In the cool CSE, several chemical processes (see Fig. 1) can alter the molecular 
abundance fractions, originally established in the dense stellar atmosphere (as a 
result of the nucleosynthesis and DU events). Hence, the absolute abundance of 
most molecules can not be used as direct tracer of nucleosynthesis and subsequent 
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DUs^. Isotopic ratios are by far the most diagnostic tracers of the stellar origin of 
elements, as they are very sensitive to the precise conditions (especially the tem- 
perature) in the nuclear burning regions, and different isotopes of the same element 
can form in very different reaction chains (e.g. Karakas et al. 2008). Moreover, 
comparisons of isotopologue fractions are less affected by systematic effects than 
are comparisons of the abundances of different molecules. There is, however, 1 
chemical process which might affect some isotopologue species differently (as CO 
and CS): chemical fractionation (Mamon et al. 1988). CO chemical fractionation 
could enhance the abundance of ^^CO in the outer envelope. However, chemical 
fractionation competes with selective photodissociation, a process that preferen- 
tially destroys ^^CO but does not significantly affect optically-thick ^^CO (Mamon 
et al. 1988). These two mechanisms are thought to usually compensate each other 
(Milam et al. 2009). 

However, observed isotope fractions for AGB and supergiant winds are (still) 
very scarce. Although, equipped with a relative small mirror compared to some 
ground-based facilities, Herschel is sensitive enough to observe some of the less 
abundant isotopologue species. Line emission of ortho- and para-water in the 
different isotopologue states H^^O, H^^O, and H^^O, as well as ^^C^^O, "C^^O, 
i^Ci^O, i^ci^O, ^^Si^^O, 29Sii6o, 30Sii^O, H^^Cl, tf'Cl, 2^Si32S, 29Si32s, 30si32s, 
28Si"S, ^SSi^^S, 13C32S, 12C33S, i^c^^s, Hi^Ci^N, Hi^Ci^N etc. are de- 

tected in the HIFI, PACS and SPIRE spectra (Royer et al. 2010; Decin et al. 
2010d,b; Schoier et al. 201 1 ; Agundez et al. 201 1). Often, one can only deduce a 
lower limit to the isotopologue fraction^ due to opacity effects. Analysing several 
lines covering different excitation regimes partly solves this issue. 

For the red supergiant VY CMa Royer et al. (2010) derived an ortho-to-para 
ratio (OPR) for H2O of 1.27, significantly lower than the expected value of 3^^\ If 
opacity effects are not compromising these results, this reflects a spin temperature 
~15 K indicating formation of H2O in non-LTE conditions. Schoier et al. (201 1) 
derived for the S-type AGB star;^^ Cyg an ortho-to-para OPR H2O ratio of 1.4 ± 
0.34, a i^CO/i^CO ratio of 43 + 6, a Hi^CN/Hi^CN ratio of 56 + 28, a ^^SiO/^^SiO 
ratio of 8 + 2, and a ^^SiO/^^SiO ratio of 19 + 8. The quite high i^C/^^C ratio 
suggests that this ratio is higher in S-type than in M-type AGB stars, and more in 



^Some indirect hints can, however, be derived. E.g., the low HF-abundance derived for 
IRC 10216 argues towards the fact that '^F nucleosynthesis is not effective in AGB stars (Agundez 
etal. 2011). 

^defined as the fraction of the most abundant species to the least abundance species 
^reflecting a high-temperature (~50) equihbrium water formation 



21 



line with those derived for carbon stars. The Si isotopic ratios can be compared 
with the ^^Si/^^Si and ^^Si/^'^Si ratios inferred for the Sun (19 and 29, respectively), 
the M-type AGB star IK Tau (13 and 40, respectively, Decin et al. 2010d), and 
IRC +10216 (15 and 20, respectively, Kahane et al. 2000). Again, these isotopic 
ratios ofx Cyg are more in line with the carbon-rich AGB star IRC +10216. 

For another AGB target, the O-rich AGB star IK Tau, Decin et al. (2010d) 
derived in OPR for H2O of 3, and Ri^O/Hi^O = 600 + 150, H^^O/Hf O = 200 + 50 
(well below the solar values of '''0/"0~2632 and '^0/^*^0-499). The lower than 
solar '^0/'*^0 ratios cannot be explained by any stellar evolution model in the 
literature. However, IK Tau is not the only Galactic star with a low ^^O/^^O ratio: 
some barium stars analyzed by Harris et al. (1985) also have a low ^^O/^^O value. 
The observation of the higher excitation ^^SiO J=14-13 line with HIFI, resulted 
in a reduction of the ^*^SiO abundance in the inner wind with a factor 2 compared 
to the value derived by Decin et al. (2010c), proving the strength of observing 
higher-excitation lines. 

It is anticipated that the Herschel observations of other evolved stars, com- 
plemented with other ground-based data sets, will add new information to this 
discussion. If not too much plagued by opacity effects, these will yield useful 
constraints for stellar evolution models. 

3.2. Non-equilibrium chemistry in the inner envelope 

For a long time, the gas chemical composition in the CSE was believed to 
be dominated entirely by the C/0 ratio of the photosphere. A C/0 ratio greater 
than one implied that all the oxygen was tied in CO, leading to an oxygen-free 
chemistry (the so-called carbon stars), whereas a C/O ratio of less than one meant 
that no carbon bearing molecules apart from CO could ever form in an oxygen- 
rich (O-rich) environment (M-type stars). Stars having a C/O-ratio ~ 1 are called 
S-type AGB stars. This picture, based essentially on thermal equilibrium consid- 
erations (Tsuji 1973) applied to the gas, has been disproved by the detection of 
SiO at millimeter (mm) wavelength in carbon-rich (C-rich) AGBs (e.g. Bujarrabal 
et al. 1994). As for O-rich AGBs, CO2 infrared (IR) transition lines were detected 
in various objects with the Short-Wavelength Spectrometer (SWS) on board the 
Infrared Space Observatory (ISO) (e.g. Justtanont et al. 1998). 

Originally, it was thought that the carbon species observed in the winds of 
O-rich stars were produced via photochemical processes in the outer envelope. 
However, Duari et al. (1999) showed that pulsation-driven shock-induced non- 
equilibrium chemistry models predict the formation of large amounts of a few 
carbon species, Uke HCN, CS and CO2, in the inner envelope of O-rich AGB 
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stars: these molecules are formed in the post-shocked layers and are then ejected 

in the outer wind as 'parent' species. Later on, it was shown theoretically (Cher- 
chnefF 2006) and observationally (Decin et al. 2008) that the inner winds of AGB 
stars show a striking homogeneity, independent of their C/0 ratio and stage of 
evolution. 

The Herschel observations have confirmed this picture, showing that integrated 
line intensities ratios are often much larger than expected from thermodynamic 
equilibrium (TE) chemical production models (e.g. Royer et al. 2010; Decin et al. 
2010b; Schoier et al. 2011). The derived inner wind abundance might be orders 
of magnitude higher than the TE predictions. 

However, while the observations are showing overwhelming evidence for the 
occurrence of complex non-equilibrium chemistry, a direct confrontation with the- 
oretical model predictions is still lacking. The reason for this is that the non- 
equilibrium chemistry models as presented by, e.g., Agundez and Cernicharo 
(2006) and Cherchneff" (2006) only include a very restricted set of (circum)stellar 
parameters tuned toward predictions in agreement with the observations for one 
specific target. This was also the concept in the recent publication of Cherchnefi" 
(201 1). To explain the recent detection of warm water vapor in the inner envelope 
of the C-rich AGB star IRC 102 16 by Decin et al. (2010a), the chemical model 
was updated by including a greater completeness of the chemistry, and including 
a more accurate chemistry for Si and S. The impact of this updated chemical net- 
work was huge, with the water vapour abundance at ~3 R^ increasing with few 
orders of magnitude. However, it is difficult to translate this type of results to 
other evolved stars without proper knowledge of the sensitivity of the results to 
different physical, chemical and numerical input parameters. Indeed, water vapor 
is now detected in all carbon-rich AGB stars surveyed with Herschel, with the ex- 
ception of 1 target. An analysis of the water-rich sources shows a large variation 
in the derived water abundances (see Fig. 12). E.g., V Hya (with a mass-loss rate, 
M, of 1 X 10"^ Mcs/yr) and CIT 6 (M=4.2 x 10"^ Mo/yr) have a water peak abun- 
dance relative to H2, /H20, of 4.5 x 10"^ and 3 x 10^^, respectively. This is up to 2 
orders of magnitude larger than in the case of IRC -1-10216 (M=2.1 x 10"^ Me/yr, 
/h20=2 X 10"^), and a factor of 10-100 larger than the predictions of Agundez 
et al. (2010). Currently, we do not know the key parameter(s) determining the 
water abundance. 
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Figure 12: Relative line strength of some water vapor lines as observed with Herschel/PACS in 
the bright carbon sources IRC +10216, V Hya, CIT 6, and V Cyg. The line-to-continuum ratio is 
clearly highest for V Cyg and lowest for IRC -1-10216, and can vary with a factor of a few from 
source to source. 
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3.3. Dust mineralogy and gas- grain reactions 
3.3.1. Dust mineralogy 

The spectra of both AGB stars and supergiants are ideal tracers of the solid- 
state and molecular species that form in the cool CSE (see Fig. 1 1). Most of the as- 
tronomical solid state features are found in the near and mid-IR ranges. ISO, espe- 
cially with its SWS and LWS revolutionized our knowledge of dust and ice around 
stars (Waters 2004). In the LWS range, overlapping with Herschel PACS, most of 
ISO's spectroscopic dust observations were suffering from signal-to-noise prob- 
lems in all but the brightest AGB stars. The higher sensitivity of Herschel could be 
crucial in identifying/analyzing few of these solid-state features. Dust-species like 
crystalline water-ice at 61 fxra, forsterite at 69 jira, hibonite at 78 jira, and calcite at 
92.6 lira are expected or have already been detected by ISO (Molster et al. 1999, 
2002; Kemper et al. 2002; Mutschke et al. 2002). Other measured features lack an 
identification, e.g., the 62-63 jxm feature with candidate substances like dolomite, 
ankerite, or diopside (Koike et al. 2000; Kemper et al. 2002). Polycyclic aromatic 
hydrocarbon (PAH) features have been predicted to occur at far-IR wavelengths 
(Joblin et al. 2002). The detection of crystalline silicates and water ice features 
in the PACS range is an indicator of the formation of very large oxygen-rich dust 
grains in the cool circumstellar shells around some intermediate-mass transition 
sources evolving from the AGB to the PN stage. Alternatively, these oxygen-rich 
dust grains surrounded by water ice mantles may form in the coolest regions of 
long-lived circumstellar disks around (binary?) low-mass stars, which can then 
be preserved from destruction by the increasing UV field generated by the central 
star as it evolves towards the PN stage. 

A solid-state band clearly detected with PACS is the 69 jum feature of crys- 
talline olivine (Mg,Fe)2Si04. The peak and shape of the feature can be used as 
a dust thermometer, but also to assess the exact composition of the olivine grains 
(see Fig. 13 and Fig. 14). Combining the 69 //m band with the emission peaks at 
11.3 and 33.6 /im both eff'ects can be disentangled (de Vries et al. 201 1). 

The great advances hoped to be achieved into this field thanks to the higher sen- 
sitivity (and spectral resolution) of Herschel PACS/SPIRE compared to ISO LWS 
are not yet carried through. There are mainly three reasons for this. (1) A firm de- 
tection of a (new) solid-state species and a comprehensive analysis of its strength 
are closely linked to a highly accurate relative flux calibration. Early 2012, the 
relative calibration accuracy of the PACS spectrometer is typically between 5- 
10% (depending on the pointing and possible drifts), which might improve to 1% 
when the signal of diff'erent spatial pixels can be combined. (2) Thanks to the 
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Figure 13: Left-hand figure: The 69 yum band of forsterite (Mg2Si04) is of special interest because 
its opacity is sensitive to the temperature of the grains. If the grains become hotter, the band 
shifts to the red, becomes broader and it drops in strength (Suto et al. 2006). This is illustrated in 
case of a temperature of 50 (red) and 300 K (blue). Right-hand figure: Forsterite is a crystalline 
olivine ((Mg, Fe)2Si04) with a Mg-Fe composition that is purely magnesium rich. The opacity 
is dependent on the iron content (Koike et al. 2006): starting with pure forsterite and increasing 
the iron content of the crystalline olivine, the 69 //m band shifts to the red and becomes slighty 
broader. 
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Figure 14: Fits to the 69 jim band of forsterite using temperature dependent opacities (de Vries 
et al. 201 1). Fits are shown for the post-AGB star HD 161796 and the planetai-y nebula NGC 6543. 
The forsterite dust temperature is lower in HD 161796 than in NGC 6543. 
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higher spectral resolution of Herschel, it is clear that (molecular) gas lines are 
contributing to the flux emission at the same wavelengths (see Fig. 15). A dedi- 
cated analysis of these gas features is needed to properly determine the strength 
of the dust emission bands. (3) Most of the solid-state features have the largest 
extinction coefficients in the near to mid-IR wavelength range. 
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Figure 15: Comparison between the 69jum forsterite feature in the oxygen-rich AGB star 
OH127. 8-1-0.0 (M = 5 X 10"^ Mo/yr, Lombaert et al. 2012) and the rich molecular Une spectrum 
of the oxygen-rich AGB star R Dor {M - 1.2 x lO^^Mg/yr). The upper panel shows the spectra 
of both targets between 65 and 72 /im at a spectral resolution of ~1750. The lower left spectrum 
zooms into the 69 /im region. The lower right figure shows the spectrum of R Dor around 69 
micron in B2A (2nd order, black) at a resolution of ~1750 and in band 3 A (3rd order, blue) at a 
resolution of ~4500. The strongest line peaks in the band 3A spectrum between 68.5-69.25 ;um 
are due to ortho- and para-HaO, '^CO, ^'^SiO, and SO2. 
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3.3.2. Gas-grain reactions 

The abundance fractions of quite some molecules show a depletion pattern in 
the dust formation region due to their role in the formation of solid-state species. 
It is thought that some molecules, as CO, CS and HCN, are quite stable and travel 
the entire envelope unaltered until they reach the photo-dissocation region of the 
outer wind, because these molecules do not participate in the formation of dust 
grains like silicates and corundum (Duari et al. 1999). 

Thanks to the long wavelength range covered by the three Herschel spectrome- 
ters different excitation levels of several molecules can be sampled (see, e.g. Decin 
et al. 2010b) with the potential to derive possible (molecular) depletion in the en- 
velope. An example is shown for the 0-rich AGB star IK Tau in Fig. 16. SiO is 
clearly depleted in the intermediate wind region, possibly due to the formation of 
SiOi (via reaction with OH) whose condensation product is silica and which is 
also thought to play an active role in the formation of different types of silicates. 




Radius [R.] 

Figure 16: Fractional abundance stratifications for '^CO, '^CO, ^'^SiO, ^'SiO, HCN, SO, ortho- 
H^^O, para-H^^O, ortho-H^^O, para-H^^O, ortho-H^'^O, and para-H^**0 as derived by Decin et al. 
(2010d). For all molecules (except water), the full line represents the results obtained by Decin 
et al. (2010c). For '^CO and ^^SiO, the new results based on the HIFI data are shown as dotted 
lines. The fractional abundances for all water isotopologs and isomers are based on HlFl data 
alone. 
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Herschel has already observed tens of AGB and supergiant stars. However, 
only for 3 of them (IRC +10216, x Cyg and IK Tau) a detailed study of the 
molecular depletion pattern is already performed. This kind of analyses are com- 
putationally very expensive and ask quite some dedicated effort of the modeler 
to analyze tens of lines of different molecules consistently. But even then, we 
have to realize that this type of results are quite often obtained assuming a spher- 
ically symmetric geometry. The Herschel images (see Sect. 2), complemented 
with other optical and sub-millimeter data, show that this is often a too simplistic 
approach. 

Another illustrative example of the close interaction between dust and gas is 
described by Lombaert et al. (2012). In this paper, it is demonstrated how water 
vapor lines should be modeled with care, due to a degeneracy between the water 
abundance and dust-to-gas ratio (see Fig. 17). Water vapor is partly coUisionally 
excited, but for many transitions radiative excitation by the near-IR stellar con- 
tinuum photons or the thermal dust radiation field are important in establishing 
the level populations. Extra constraints on the water vapor abundance in the in- 
ner envelope can be obtained from the water ice feature at 3.1 jxm, since the H2O 
ice fraction leads to a minimum required H2O vapor abundance (Lombaert et al. 
2012). 




Figure 17: HerschelfPACS observations of OH127. 8+0.0 (black) compared to model predictions. 
Both a model with a low dust-to-gas ratio if/ of 0.003 and a high water abundance of 1 x 10""* (blue 
full line), as well as a model with a high dust-to-gas ratio of 0.01 and a lower water abundance of 
5 X 10"^ (red dashed line) can predict the PACS H2O lines equally well. 
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3.4. Photodissociation of molecules and ion-molecule reactions 

The penetration of highly energetic interstellar UV and cosmic ray radiation 
in the tenuous outer envelope induces the photodissociation of so-called 'parent' 
molecules^, which in its turn leads to the formation of new 'daughter' species. 
E.g., in the outer envelope of both C- and 0-rich sources, CN is thought to be 
produced by the photodissociation of HCN and neutral-neutral reactions (Willacy 
and Millar 1997; Decin et al. 2010c). 

The new data obtained with Herschel show indications that this process of pho- 
todissociation and the production of new molecules due to ion-molecule reactions 
might even be an important chemical player in the inner envelope. One key dis- 
covery made by Herschel is the detection of warm water vapor in the inner enve- 
lope of the carbon-rich AGB star IRC -1-10216 (Decin et al. 2010a; Neufeld et al. 
2011b). The existence of water vapor around carbon-rich AGB stars is remark- 
able, because, unlike those of oxygen-rich stars, the photospheres of carbon-rich 
stars are expected to contain very little H2O as they are dominated by CO, HCN, 
and C2H2 under conditions of thermochemical equilibrium. The water abundances 
derived from Herschel observations of carbon-rich AGB stars are typically 3 to 4 
orders of magnitude larger than the expected photospheric abundance. At least 
in the case of IRC -1-10216, Herschel's discovery of water in the warm inner en- 
velope imposes significant observational constraints on the various theories that 
have been proposed. Thus far, at least two possible origins appear to be consistent 
with the available data. In one picture (Decin et al. 2010a; Agundez et al. 2010), 
the water is released by the photodissociation of ^^CO and SiO, thanks to the pen- 
etration of interstellar UV radiation through the clumpy outflow surrounding the 
star. In another scenario (Cherchnetf 201 1, see Sect. 3.2), water results from non- 
equilibrium chemistry associated with pulsationally-driven shock waves. Other 
origins including (1) the vaporisation of icy objects (comets or dwarf planets) in 
orbit around the star (Melnick et al. 2001); (2) Fischer- Tropsch catalysis on the 
surface of small grains (Willacy 2004); or (3) photochemistry within an outer, 
photodissociated shell (Agundez and Cernicharo 2006) could be ruled out since 
they predict an absence of abundant water within ~100 AU of the star, which is 
not compliant with the relative strength of the observed higher-excitation water 
vapor transitions. 

A second important detection has been the widespread presence of water va- 



i.e. formed in the outer stellar atmosphere or inner envelope and ejected into the intermediate 
envelope 
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por in carbon- rich AGB stars. A small Herschel/HLFl survey for water vapor in 

eight carbon-rich AGB stars obtained a water detection in every star surveyed, 
suggesting that IRC +10216 is unusual only in its proximity to Earth, and that 
the presence of water in carbon-rich AGB stars is nearly universal (Neufeld et al. 
201 la). In 6 (out of 7) sources already observed with PACS, several low, mid, and 
sometimes higher excitation H2O lines were detected (see, e.g.. Fig. 12). How- 
ever, much to our surprise, no water lines were detected in 1 target. While the 
results of Herschel obtained in its first year of observations surprised us in terms 
of the detection of warm water vapor in carbon-rich stars, the observations taken 
during the next year implied an inversion of the science question, i.e. 'Why don't 
we see water in all carbon-rich sources?'. The key parameter determining the 
water abundance is currently not known. 

3.5. CSM-ISM interaction region 

The detection by Herschel of the ubiquitous nature of a (turbulent) astropause 
at the outer edge of the circumstellar envelope stimulates the wish to unravel the 
chemical composition in the CSM-ISM interaction region. In the end, the ther- 
modynamical and chemical structure in this zone determines the chemical yields 
injected into the ISM. Currently, the knowledge on the chemical content is very 
limited. The bow shock regions are either detected in the infrared or far-UV (e.g. 
Martin et al. 2007; Ueta et al. 2006; Cox et al. 2012). 

The origin of the FUV/NUV emission of the bow shock around o Cet has 
been attributed to H2 molecules in the cold gas which are coUisionally excited 
by hot electron from the postshock gas (Martin et al. 2007). The origin of the 
FUV/NUV emission in the interaction zone between the ISM and the stellar wind 
of CW Leo has not been modeled in detail, although it was suggested by Sahai 
and Chronopoulos (2010) that also here coUisionally excited H2 emission might 
be the cause. The FUV/NUV emission around CW Leo peaks at slightly larger 
distances from the central target than the PACS/SPIRE infrared flux excess (Lad- 
jal et al. 2010). Note that not all sources showing a clear IR bow shock detection 
exhibit UV emission. An illustrative example in this case is a Ori (Decin et al. 
2012a). 

The infrared emission is probably due to thermal emission of cold dust compo- 
nents in the interaction region, with temperatures around 30-150 K. An attempt to 
detect probable contribution from low-excitation atomic lines, such as [01] 63 yum 
and [Cn] 158 //m, at the IR wavelength bands, has been unsuccessful in case of 
a Ori (Decin et al. 2012b). Neither is the bow shock detected in low-excitation 
CO rotational lines. Gas temperatures calculated by hydrodynamical models (see 
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Sect. 2.3) are in the order of few thousand K in the shocked regions. Currently, 

it is absolutely unclear how this high gas temperature will impact the atomic and 
molecular abundance pattern in this region. In addition, we still do not know how 
the turbulence in the astropause will impact the grain temperature. 

4. Some old and new questions 

Since the first detection of a thermally excited circumstellar line, CO v = 0, 7 = 
1 - at 2.6 mm, toward IRC 10216 by Solomon et al. (1971) our knowledge on 
the cool envelopes surrounding late-type LIMS has increased steadily. The detec- 
tion in 1971 was interpreted as being a measure of the radial velocity of a steady 
mass-loss process, or alternatively, as being the result of an explosive event which 
ejected a shell from the central star. Later on, more and more observations cor- 
roborated the interpretation of an important mass-loss process during the AGB 
and supergiant phase. The idea of a steady mass-loss process has remained for a 
long time, and often a constant mass-loss rate was assumed in the modeling for 
the sake of simplicity. That way, the circumstellar envelopes could be modeled 
using ID codes. However, astonishing evidence is again given by the Herschel 
observations that these envelopes are far from smooth. Complexity is the rule, not 
the exception. Since these complex structures are the result of a mass-loss process 
initiated in the upper stellar atmospheres, the upper atmosphere structure comes 
into question. For supergiants, the huge convective cells detected in the upper 
atmosphere (Lim et al. 1998; Chiavassa et al. 2010) are indeed the direct proofs 
of non-homogeneous atmospheric structures. Cool gas should co-exist with the 
hotter gas. Related to this convective motion, is the idea of a solar-like magnetic 
activity cycle (Soker 2000). Owing to the appearance of magnetic cool spots dur- 
ing the active phase, the gas temperature can be locally reduced. This can enhance 
dust formation and thus lead to a higher mass-loss rate from magnetic cool spots. 
For the less massive counterparts, the AGB stars, radiation pressure on dust grains 
is still thought to be the main driver of the mass-loss process. However, in that 
case, one would expect a more or less smooth envelope structure. The complex 
inner envelope structures shown in Fig. 3 prove that another physical process, of 
which the strength varies across the upper stellar atmosphere and/or inner enve- 
lope, can not be neglected. A possible suggestion is related to dust shadowing 
or multi-scattering processes, but the data currently lack spatial resolution to cor- 
roborate this idea. How this affects the wind acceleration should still be sorted 
out. 
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We also should realize that stellar structure and evolution models are still us- 
ing a simplified (analytical) description of the mass-loss rate. The Herschel ob- 
servations show that one might start wondering about this type of assumptions. 
Improving this description is important since the mass-loss rate determines the 
fraction of stars which finally will end their Uve as supemovae or white dwarfs. 

In my opinion, one of the main answers we can get from the Herschel obser- 
vations concerns the empirical determination of the global envelope structure, i.e. 
temperature T{r) and abundance structure e(r). Thanks to the extended wave- 
length coverage of the three instruments, we are sampling the CO ladder from 
/ = 4 - 3 to / = 53 - 52, even in some rare cases seeing rotational transitions in 
vibrationally excited states. Complemented with some ground-based data (to get 
the lower rotational transitions), these CO data form an excellent thermometer. 
Last few years, effort has been spent to derive the temperature structure from first 
principles (see, e.g., Decin et al. 2006). However, still quite some assumptions are 
part of the functions determining the cooHng and heating rates per unit volume. A 
more reliable temperature structure can be derived from the data themselves. The 
same holds for the abundance structure. The Herschel spectra are extremely rich, 
especially in terms of molecular line transitions. In the case that different excita- 
tion regimes are sampled 6(r) can be derived, yielding observational constraints 
for chemical models. For the first time, infrared and submillimeter data of a large 
set of evolved LIMS, covering different (circum)stellar parameters, is available to 
the community. A homogeneous analysis of these data will form an important 
testbed for the chemical models. The observation of water vapor in carbon-rich 
stars can serve as a good example. While in 2001, it was argued that water vapor 
in this type of carbon-rich environment was due to the vaporization of icy bod- 
ies (Melnick et al. 2001), the detection of high-excitation water vapor Unes by 
Herschel gave arguments against this interpretation (Decin et al. 2010a). At that 
moment, the idea of photo-induced chemistry in a non-homogeneous environment 
was postulated. Later on, updated chemistry models showed that pulsationally- 
induced non- equilibrium chemistry could explain the observations as well (Cher- 
chneff 2011). And while in 2010, we could prove the ubiquitous nature of water 
vapor in C-rich envelopes, in 2012 we start wondering why we do not see water 
vapour in all C-rich AGB stars (see Sect. 3.4). This proves that in the last ~10 
years time our ideas on the chemistry in these envelopes have changed drastically. 

A new field nicely uncovered by Herschel concerns the CSM-ISM interaction 
region. The data shown by Cox et al. (2012) show for the first time a whole zoo 
of interaction phases, roughly sub-divided in four categories. The complexity at 
the contact discontinuity is huge, the dynamical and temperature structure still 
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largely unknown, let alone for the chemical pattern. How does this harsh en- 
vironment influence the structure of the grains and the photodissociation of the 
molecules? And how does this energetic environment determine the composition 
of the chemical yields enriching the ISM? 

5. Future perspectives 

In my opinion, two keywords on the exploitation of the Herschel in the next 
10 years are 'confrontation' and 'consistency', if possible coupled to each other. 
We are still in the first exploration status of the Herschel data, scrolling through 
images and spectra and looking for spectacular results. These first results are im- 
portant in their own right, and do give us already some nice clues on late stages 
of stellar evolution. A next step should aim for an in-depth confrontation of the 
data with current theoretical models. On the one hand, this will demonstrate the 
strength of the models, and show how some general quantities might be extrap- 
olated to other evolved stars, maybe in other galaxies. On the other hand, this 
will inevitably show us the limitations of the models/our knowledge. One of the 
main limitations the field of evolved- stars research is currently confronted with 
concerns 'consistency' between applied models. There are models focusing on 
the gas chemistry or on grain chemistry, models solving the thermodynamic equa- 
tions, models focusing on multi-D radiative transfer etc. However, none of the 
current models is fully consistent. The dust grains and gas molecules are intri- 
cately interacting in the CSEs, both radiatively and dynamically. Therefore it is of 
interest to adopt a model that is as self-consistent as possible. One aspect concerns 
the fact that AGB stars are dynamically unstable, resulting in the star changing its 
luminosity when going from maximum to minimum light. This severely aff'ects 
its spectral energy distribution, as shown by e.g. De Beck et al. (2012). However, 
not only the thermal dust emission is affected, but a change in luminosity (and 
hence dust opacity) will change the line strengths of quite some lines, both in 
the infrared (see, e.g., van Malderen 2003) and sub-millimeter regime (see, e.g., 
Lombaert et al. 2012, and Fig. 18). This eff'ect should be taken into account when 
modeling data taken at different phases. 

Another example concerns the formation of some gas species (as H2O), which 
might influence the thermal structure (cooling or heating the envelope). This local 
temperature change might afi'ect in its turn the formation of specific dust grains, 
which are known to act as wind drivers. Added to this are effects as backwarming 
or dust shadowing in combination with, e.g., multi-scattering of stellar and dust- 
emitted photons. This might yield non-isotropic dust formation. Admittedly, this 



34 




Figure 18: Comparison between theoretical spectra for different CO transitions for a model with 
parameters similar to the 0-rich giant IK Tau (Decin et al. 2010c). The full black line represents a 
model for a luminosity of 1000Lo(TefF=2200K, R+=1.5x lO'^ cm), the red dotted line is the result 
for a model with a luminosity of 2000 Lq, simulated by increasing the effective temperature, Tefj, 
to 2600 K, and lowering the stellar radius to R,^=1.3 x lO'"* cm. A change in luminosity mainly 
affects the higher-excitation lines formed closer to the central star. 
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type of 'consistent' models are CPU consuming, but the theoretical knowledge 

on each specific sub-field is there. It just should be combined. While Herschel 
gives us a whole zoo of infrared and submillimeter data, the theoretical models 
are currently lagging behind. 

Another issue, with which the Herschel data are confronting us, is our limited 
knowledge on the chemical properties of many gas and dust species. Even in case 
of relatively 'simple' molecules, as C2H our radiative transfer treatment is still 
very limited: C2H is a linear molecule with an electronic ground-state config- 
uration of which nowadays radiative ransfer calculations are done assuming a 
approximation. The dipole moments of the fundamental vibrational states of C2H 
are still unknown (De Beck et al. 2012). This lack of knowledge might severely 
limit the accuracy of the derived results. This has been demonstrated in the case 
of H2O, of which it is now well established that one overestimates the H2O abun- 
dance by an order of magnitude in the case one neglects the infrared pumping of 
H2O via various vibrational bands (Agundez and Cemicharo 2006). This situa- 
tion will hopefully improve in the near future for simple 2 and 3 -atomic molecules 
thanks to initiatives as the Cologne Database for Molecular Spectroscopy (CDMS, 
Miiller et al. 2005)'', the JPL catalog (Pickett et al. 1998)^°, basecol (Dubemet 
etal. 2006) etc. 

6. Conclusions 

To conclude, the harvest by the Herschel observations in the field of evolved 
stars is large. New ideas on the chemical, thermal and dynamical structure in the 
envelopes created by the dominant mass-loss process have emerged. As of today, 
this is reflected in more than 150 refereed-joumal papers in the field of evolved 
AGB and supergiants stars containing Herschel data. Quite some evolved stars 
are scheduled to be observed with Herschel. The comprehensive analyses and 
consistent confrontations with theoretical models will still take years. 

A next observational step, which has already begun, concerns high spatial res- 
olution observations, as will be provided by instruments as ALMA. With its high 
sensitivity, ALMA will enable us to resolve the inner wind structures of the most 
nearby targets and the high-resolution maps of the winds will be an excellent tracer 



'http://www.astro.uni-koeln.de/cdms/ 
^''http://spec.jpl.nasa.gov 
^ ^http://basecol.obspm.fr/ 
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of the dust and gas distributions, and will in that way complement the extremely 
valuable Herschel data. 
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